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ABSTRACT 

Using star-forming galaxies sample in the nearby Universe (0.02<z<0.10) selected from the 
SDSS (DR7) and GALEX all-sky survey (GR5), we present a new empirical calibration for 
predicting dust extinction of galaxies from Ha-to-FUV flux ratio. We And that the Ha dust 
extinction (Anct) derived with Ha/H/3 ratio (Balmer decrement) increases with increasing 
Ha/UV ratio as expected, but there remains a considerable scatter around the relation, which 
is largely dependent on stellar mass and/or Ha equivalent width (EWhq). At fixed Ha/UV 
ratio, galaxies with higher stellar mass (or galaxies with lower EWna) tend to be more highly 
obscured by dust. We quantify this trend and establish an empirical calibration for predicting 
Ana with a combination of Ha/UV ratio, stellar mass and EWho, with which we can suc¬ 
cessfully reduce the systematic uncertainties accompanying the simple Ha/UV approach by 
^15-30%. The new recipes proposed in this study will provide a convenient tool for predict¬ 
ing dust extinction level of galaxies particularly when Balmer decrement is not available. By 
comparing Ana (derived with Balmer decrement) and Auv (derived with IR/UV luminosity 
ratio) for a subsample of galaxies for which AKARIEIR photometry is available, we demon¬ 
strate that more massive galaxies tend to have higher extra extinction towards the nebular 
regions compared to the stellar continuum light. Considering recent studies reporting smaller 
extra extinction towards nebular regions for high-redshift galaxies, we argue that the dust 
geometry within high-redshift galaxies resemble more like low-mass galaxies in the nearby 
Universe. 
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1 INTRODUCTION 

Star formation rate (SFR) is one of the most fundamental param¬ 
eters which characterise the nature of galaxies. Deriving SFRs of 
galaxies is therefore an important task in the extra-galactic astron¬ 
omy. There are various indicators proposed for star-formation ac¬ 
tivity in galaxies, including rest-frame UV continuum light, nebu¬ 
lar emission lines (such as Ha), mid- to far-infrared du st thermal 
emission, and radio continuum l uminosity (see review bv iKennicutll 
ll998l : lKennicutt & Evanslbol^ . 

When measuring SFRs, dust extinction correction is always 
an important issue. Particularly in the case of star-burst galaxies, it 
is possible that only a tiny fraction of UV light can escape from the 
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galaxy because they are heavily obscured by dust. This leads to an 
order of magnitude extinction correction for estimating the intrinsic 
UV luminosity (hence SFRs). For studies of high-redshift galaxies, 
where the cosmic star formation rate density is an order of magni¬ 
tude higher than the prese nt-day universe fe.g. lMadau et alj|l996l : 
iHopkins & Beacoral200d) . luminous infrared galaxies (LIRGs) or 
ultra-luminous infrar ed galaxies (UUIRGs) are much more com - 
mon population (e.g. lUe Floc’h et aflFoOSl : iMagnelli et alJl2oTTh . 
and therefore the effect of dust must be carefully taken into ac¬ 
count. Far-infrared (FIR) observations, on the other hand, allow a 
direct measurement of dust thermal emission, and in particular, ob¬ 
servations covering a peak of spectral energy distribution (SED) 
of galaxies (usually located at around Arest~100 /im) allow us to 
derive total IR luminosity (hence dust-enshrouded SFR) with rea¬ 
sonable accuracy. A problem is that such FIR observation requires 
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space telescopes. Unfortunately, due to the poor spatial resolution 
of the IR space telescopes ever launched, the depths of FIR observa¬ 
tions are always limited by source confusion, making it impossible 
to detect individual galaxies at high-redshift, except for exception¬ 
ally luminous objects. 

Another important tool for measuring SFR is HaA6563 line, 
which is well-calibrated in the local Universe. The Ha emission is 
emitted in star-forming Hit regions (near short-lived, O-type stars), 
and so the Ha luminosity is expected to be proportional to on¬ 
going SFR. The Ha line is located at optical wavelength (hence 
less sensitive to dust extinction effects than UV light), and there¬ 
fore Ha line is recognised as an excellent indicator of SF activity. 
Of course, Ha line is also affected by dust. Even for local spiral 
galaxies with moderate levels of star formation , the dust extinc - 
tion at Ha is not negligible (Ahu ~1 mag: e.g. lKennicutill983h . 
In the case of more active starbursts, the Ha line is reported to 
be more heavil y obscured by dust, an d Aua exceeds ~3-m ag in 
extreme cases l Poggianti & WuIItOOOI : iKovama et alj|2010h . It is 
therefore ideally required to combine the H/3A4861 line flux, with 
which one can derive dust extinction level using Ha/H(l ratio (i.e. 
Balmer decrement). This method needs to assume the shape of ex¬ 
tinction curve, as well as the electron density and temperature, but 
other than that, it allows us to derive dust extinction levels purely 
based on physics. With recent advents of sensitive NIR instruments 
installed on 8-m class telescopes, it is becoming easier to obtain 
Ha information for high-redshift (z > 1) galaxies; e.g. with spec¬ 
troscopy and/or narrow-band imaging. However, H(1 line is usu¬ 
ally much fainter than Ha, making it very challenging to apply the 
Balmer decrement method for measuring dust extinction for high- 
2 galaxies. Indeed, it is becoming possible to obtain unprecedent¬ 
edly high-quality NIR sp ectra for high-z galaxies (e.g. lSteidel et al.l 
l2014l ; lReddv et alj20la) . but even with such high-quality data, it is 
still challenging to detect H/1 lines from individual galaxies. 

Our idea is to use observed Ha-to-UV flux ratio to predict 
dust extinction levels of galaxies. For high-redshift studies, rest- 
frame UV flux density can easily be traced by optical photom¬ 
etry. Considering the fact that Ha is less sensitive to dust ex¬ 
tinction than UV continuum light, it is expected that the Ha/UV 
flux ratio can provide a crude test for dust extinction level (i.e. 
galaxies with higher Ha/UV ratio are expected to be dustier). 
It is of course true that an intrinsic Lua/Luv ratio can change 
with galactic age, SF history, metallicity, or initial mass function 
(IMF), and therefore it is not straightforward to directly link the 


Ha/UV r a tio to th e dust extinction properties dWuvts et al.l 


Lee_et_aJ 2009 ; |Pflamm-Altenburg, Weidner. & Krour ^ 


2013; 


200S; 


Zeimann et alj|2014l) . Nevertheless, as demonstrated by lBuat et al.l 


l 2002h with a small sample of local star-forming galaxies, the 
Ha/UV ratio is indeed well correlated with dust extinction proper¬ 
ties. In this paper, we will revisit this issue and attempt to construct 
an empirical calibration to predict dust extinction with Ha/UV ra¬ 
tio, by compiling a statistical sample of nearby galaxies (z < 0.1) 
drawn from SDSS, GALEX, and AKARI. 

This paper is organised as follows. In Section 2, we present 
our datasets and summarize physical quantities used in the paper. 
We use SDSS, GALEX, and AKARI data to obtain UV-, Ha-, and 
FIR-based SFRs. Our main results are shown in Section 3. We 
first present a positive correlation between dust extinction (Ahq) 
and Ha/UV flux ratio (Section 3.1), and then we draw the scatter 
around the Ana-Ha/UV correlation as functions of stellar mass 
and/or Ha equivalent width (Section 3.2-3.4). We establish an use¬ 
ful prescription for predicting dust extinction in the absence of H/3 
line or FIR photometry (Section 3.5). In Section 4, we present re¬ 



Figure 1. Summary of our sample selection and the number of the galaxy 
sample at each step during the cross-identification between SDSS, GALEX, 
and AKARI sources. 


lated analyses and discuss the dust properties of star-forming galax¬ 
ies. We first provide a careful investigation of the systematic un¬ 
certainties associated with aperture correction for SDSS data (Sec¬ 
tion 4.1). We then discuss if we can really derive the dust extinction 
levels of galaxies by using stellar mass alone (Section 4.2). We also 
discuss the relation between dust extinction and metallicity in Sec¬ 
tion 4.3. Furthermore, in Section 4.4, we discuss the extra extinc¬ 
tion towards nebular regions, and report its dependence on stellar 
mass and EWhq. Finally, we will discuss if our new recipes can be 
applicable to high -2 galaxies using our high -2 galaxy sample (Sec¬ 
tion 4.5). Our conclusion is given in Section 5. Throughout the pa¬ 
per, we adopt Dm = 0.3, Da = 0.7, and Hq = 70 km s“'^Mpc~^. 


2 DATA AND SAMPLE SELECTION 
2.1 SDSS data 

We use the spectroscopic catal ogue of the Sloan Digi tal Sky Survey 
(SDSS) Data Release 7 (DR7; [^azaiian et al.ll200^ . The spectro¬ 
scopic measurements are performed by the Max Planck Institute 
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Figure 2. Difference between total (petrosian) magnitudes and SDSS fibre magnitudes at r-band as functions of stellar mass, SFRua, EWhu. and Aho (from 
left to right). There seems to be a decreasing trend of the magnitude difference with EWjjai which will be explained by the fact that high-EW galaxies tend 
to be more compact, whilst the correlation between the magnitude difference and other galaxy properties are much weaker. To create this plot, we divide each 
panel into 90x90 sub-grid, and count the number of galaxies in each pixel. Note that we basically follow this strategy when we show the grey-scale or colour 
image plot in the remainder of this paper, unless otherwise stated. 


for Astrophysics and Johns Hopkins University (MPA/JHU group), 
and we make use of the “value-added” catalogues retrieved from 
their public websit^B The catalogue contains a total of 927,552 ob¬ 
jects, of which 868,858 are identified as unique sources: we here 
perform internal matching using 2" search radius to avoid dupli¬ 
cated objects. Furthermore, we restrict the sample with a redshift 
range of 0.02<2:<0.10, and exclude those having uncertain red- 
shifts (by applying Z_WARNING=0). With these criteria, we have 
selected 344,425 objects. The redshift range adopted here was cho¬ 
sen so that we can study galaxies in a wide range in stellar mass, 
and at the same time we can minimise a potential effect of redshift 
evolution of galaxy properties within the sample. We note that a 
small change of the redshift range applied here does not affect our 
results. 

We select star-forming galaxies using the ’BPT’ diagram 
( [Baldwin. Phillips, & TerleviM Il98lh . The BPT diagnostics re¬ 
quires [OlIl]A5007/H/3 and [Nll]A6584/Ha line flux ratio, and 
therefore these four major emission lines need to be detected. We 
first request S/N(Ha)>10 (as Ha line is usually the strongest), 
yielding 196,073 sample. We also request the following crite¬ 
ria for S/N ratio of the other three emission lines: S/N(H/3)>3, 
S/N([Nll])>3, S/N([OlIl])>2. We note that [Olll] line is the weak¬ 
est in most cases, and so we accept the [Olll] line detection down 
to S/N>2. After applying all these criteria, 177,033 galaxies are se¬ 
lected (i.e. >90% of the Ha-selected galaxies (with S/N>10) have 
significant detection at all the other three lines). Following the rec¬ 
ommendation by the MPA/JHU group, we scaled the uncertainty 
of Ha, np, [Nil], and [Olll] by 2.473, 1.882, 2.039, and 1.566, re¬ 
spectively. We also note that the line flux measurements provided 
in the MPA/JHU catalogue are continuum-subtracted (and also cor¬ 
rected for Galactic reddening), and so the effect of stellar absorp¬ 
tion is properly taken into account. _ 

Following the prescription bv iKauffmann et aT] ( l2003h . we se¬ 
lect SF galaxies with: 

log [OHIJ/H/3 < 0.61/(log([NII]/Ha) - 0.05) + 1.3, (1) 

and we also apply [Nll]/Ha<0.6 to remove any Seyfe rt galaxies 
or LINERs dKauffrnann et al.ll2003l : lKewlev et al.ll2006h . We have 
now selected 145,923 star-forming galaxies. 


^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 


Our final requirement for the sample selection is the avail¬ 
ability of stellar mass estimates (M*). Stellar mass of the SDSS 
DR7 galaxies are computed by the MPA/JHU group, by fitting 
to the SDSS bro a d-ban d ph otometry following the philosophy of 
iKauffmann et al.l ( l2003h and lSalim et alj l l200'il) . In this work, we 
use galaxies with log(M*/MQ)>8.5. Our results do not change 
even if we do not apply the stellar mass cut, but we note that the 
sample size significantly decreases below log(M*/M0)=8.5 (by 
a factor of ~20x compared with those at the peak of the M* 
distribution), and therefore we decided not to use galaxies below 
log(M*/MQ)=8.5 for statistical analyses presented in this work. 
Overall, our final SF galaxy sample contains 135,813 galaxies (with 
the median redshift of z = 0.063). Our sample selection procedure 
described here is also outlined in Fig.[T] 

We derive Ha dust extinction (Ahq,) for each galaxy using 
Ha/H/3 flux ratio (i.e. Balmer decrement) using the following equa¬ 
tion: 


tIhck — 


2.5A:Ha 
kua — fcH/3 


log 


2.86 A 

Fua/Fap) 


( 2 ) 


where the term 2.86 is the intrinsic Ha/H/J flux ratio for Case B 
recombination at a temperature of 10^ K and at an e lectron den- 
sity of rie=10^ cm~ ^. Following iGam & Best! ( |201(]|) . we assume 
ICalzetti et alj ( l200(]h dust extinction law to calculate A:hq and fcn/s, 
so that Aua is derived in a more specific form of: 


Ahc = 6.53 log {Fhc/Fhp) - 2.98, (3) 


with the median uncertainty of ~0.24 mag (based on the un¬ 
certainties regarding the Ha/H/3 measurements). We then esti¬ 
mate star formation rate (SFR) of each galaxy based on its 
Ha luminosity. After correcting the dust extinction using the 
Alia value derived above, we comp ute SFRHa using th e stan¬ 
dard calibration of iKennicutJ ( Il99^ for iKroupal ( 1200ll) IMI0: 
SFRhq,=5.5x 10“"^^Lhq, [ergs“'^]. We chose Kroupa IMF to be 
consistent with stell ar mass estimates for SD SS galaxies performed 
by MPA/JHU group dKauffrnann et al.l2003h . but the choice of IMF 
does not affect our conclusion as our results are mostly based on the 
ratio of SFRho and SFRuv- 


^ iKennicutt et alj d2009ll note that the zero point of the SFR needs to be 
reduced by a factor of 1. 44 in the case of Kroupa I MF, compared w ith the 
original lKennicutll jl998h calibration (which assumes ISalpetejfTgBu MF). 
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Figure 3. Comparison between the observed Ho-derived SFRs (with aper¬ 
ture coiTection) and UV-based SFRs for our SDSS-GALEX sample. The 
dotted line shows the one-to-one coiTelation; SFRjja obs=SFRFuv obs- 
A systematic offset towards higher SFRjja reflects the different level of 
dust extinction for Fla and UV (i.e. H« is less sensitive to dust extinction). 


We note that SDSS spectroscopic measurements were per¬ 
formed with a limited size of fibre (3-arcsec in diameter). We there¬ 
fore apply an aperture correction based on the difference between 
total and fibre-based r-band magnitudes: i.e. we assume that the 
emission-line profile follows that of continuum light. We believe 
that this assumption is reasonable, as supported by previou s studies 
jBrinchmann et al.]|2004l : iKewlev. Jansen. & Gelletll2003) . and us¬ 
ing our own sample, we also check the difference between fibre- 
and total-magnitudes (at r-band) as functions of various galaxy 
properties (Fig. IS. We find that the fibre magnitudes sometimes 
trace only 10% of total light for a fraction of galaxies, but we also 
find that the level of aperture correction does not strongly corre¬ 
late with physical properties (M*, SFR, Aua)- An exception is that 
there is a decreasing trend of the magnitude difference with EWhh, 
which is probably because galaxies having higher EWhq tend to 
be more centrally concentrated. We also verify that our aperture- 
corrected SERhq show good agreement with those derived from 
AKARI FIR photometry (see also Section 2.3). Thus, we believe 
that the aperture correction does not make any artificial trend, but 
the results presented in this work need to be verified independently 
by using large aperture photometry for measuring Ha emis- 
siorlj. More detailed discussion regarding the aperture correction 
will be provided in Section. 4.1. 


2.2 GALEX data 

The goal of this paper is to draw dust extinction properties of 
galaxies as a function of Ha/UV ratio. We here perform cross- 

® We speculate that the typical uncertainties regarding the aperture correc¬ 
tion would be ~0.15 dex (or ~0.4 mag), by computing the standard devi¬ 
ation of the difference between the aperture correction values measured at 
u-band and r-band (i.e. Art — Ar: see also Section 4.1). However, it is not 
straightforward to determine the uncertainties associated with the aperture 
correction, and we do not consider uncertainties regarding aperture correc¬ 
tion in the remainder of this work. 


identification between our SDSS galaxy sampl e with Galaxy Evo¬ 
lution Explorer (GALEX; iMartin et al.l l2005h UV sources. We 
use unique GALEX sources selected from All-sky Imaging Sur¬ 
vey (AIS) of GALEX fifth data release (GR5), reaching down 
to the depth of ~19.9 and ~20.8 [AB mag] at far-UV (FUV; 
Aeff = 1516 A) and near-UV (NUV; Aeff=2267 A) respe ctively. 
Here we use the catalogue published bv iBianchi et al.l(l201 ll) . avail¬ 
able on the MAST web sit^H We cross-identified our SDSS (spec¬ 
troscopic) galaxy sample with GALEX sources using 3” search ra¬ 
dius. We find that 92,561 SE galaxies (out of 135,813 galaxies se¬ 
lected in Section 2.1) have UV counterparts, among which 78,731 
galaxies are detected at both FUV and NUV (because of the limited 
depth at FUV). In this paper, we use star-forming galaxies detected 
at both FUV and NUV as the “SDSS-GALEX” sample. 

We apply the Galactic extinction corrections for th e UV pho¬ 
tometry using the lSchlegel, Einkbeine r. & Davis i .l99^ dust map 
and th e Galactic extinction curve of ICardelli, Clayton. & Mathid 
d 19891) for Rv ~ 3.1: specifically, we apply Apuv = 7-9E{B — 
V) and Anuv = 8.0E{B — V). We also apply fc-correction 
for t he UV photometry using publicly availab le fc-correction 
tool jChilingarian, Melchior, & ZolotuklMl2010l : 2012). We use 
FUV—NUV and NUV—r colours to predict fc-correction of indi¬ 
vidual galaxies at FUV and NUV, respectively. In this paper, we use 
FUV photometry to compute UV-based SFRs, because FUV pho¬ 
tometry is expected to be more reliable than NUV because fluxes 
at NUV wavelength are often contributed by stars with life time 
of >100 Myr. We note that the typical fc-correction value at FUV 
band turns out to be only ~0.02-0.03 mag level, and so the effect 
on our results is negligible. _ 

We derive SFRuv following the standard iKennicutJ ( Il998h 
calibration assuming the Kroupa IMF to be consistent with the 
Ha SFRs: SFRuv=9.7x I0 "^®L^,fuv [ergs“^ Hz"i]. We stress 
again that the choice of IMF would not strongly affect our re¬ 
sults since we discuss the ratio of Ha- and UV-based SFRs (and 
the focus of this paper is mainly on massive star formation where 
most IMFs agree). In Fig.[3 we compare the observed SFRho and 
SFRfuv (i-e- without dust extinction correction for both SFRs). 
There is a general trend that SFRho tend to be higher than 
SFRfuv, which is likely because of the dust extinction effect, as 
will be discussed later in this paper. 


2.3 AKARI data 


We also match our SDSS-GAL EX sample with the AKARI EIS 
bright source catalogue (ver. 1: 1 Yamamura et all l2009l : 2010) to 
obtain FIR fluxes for individual galaxies. The AKARI satellite 
dMurakami et alj l2007h is a Japanese infrared space telescope, 
which performed all-sky su r vey in the MIR and FIR wavelengt h 
range llshihara et al.l 1201 Ol : IPoi et al] 12015l : iTakita et"^ I2015I) . 
Sources listed in the AKARI FIS bright source catalogue is se¬ 
lected and flux-limited at 90/rm. The catalogue includes photo¬ 
metric information at four FIR bands (60, 90, 140, I60/rm). Their 
5 c 7 sensitivity at each band is 2.4, 0.55, 1.4, and 6.3 Jy, for 
N60(60^m), Wide-S(90^m), Wide-LfMO^tm), and NI60(I60Atm), 
respectively. Most of the AKARI FIR sources are expected to 
be Galactic objects, but it also contains a substantial number of 


Goto et al.ll2011alb; Yuan et al.l2012: 

Toba et al.ll2013l: Totani et al.l 

2014: Kilerci Eser. Goto. & Doj2014 

). Considering the PSE size of 


http ://archive. stsci. edu/prepds/bianchi_gr5xdr7/ 
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Figure 4. SFR from UV+IR photometry plotted against Ha-based SFRs 
(with aperture and dust extinction correction) for all SF galaxies in our 
AKARI-SDSS-GALEX sample, demonstrating a good agreement between 
the two independent measurements of SFR, and supporting our procedure 
for Ho aperture and dust extinction correction works reasonably well. 

~40 arcsec of the AKARI 90/im data, we search counterparts for 
the AKARI FIS sources within our SDSS-GALEX catalogue (i.e. 
78,731 star-forming galaxies at 0.02 < 2 ; < 0.1; see Section 2.2) 
using 20” radius, and find 1200 AKARI/FIS sources have counter¬ 
parts in the SDSS-GALEX catalogue with the median redshift of 
2=0.038. In this work, we use this 1200 SE galaxies as our final 
“AKARI-SDSS-GALEX” sample. 

Because of the limited depths of AKARI all-sky survey, most 
galaxies in the AKARI-SDSS-GALEX catalogue are not detected 
at all four AKARI FIS bands. We therefore derive total infrared 
luminosity (Lir) using the WIDE-S (9Q/xm ) and WIDE-L (140pm) 
photometry as presented bv lTakeuchi et al.l ( l2010h : 

log Lir = 0.964 log ifSii + 0.814, (4) 

where LI'^kari = Ai/9o^niL./(90pm) -|- Atyi40MmL./( 140pm). 
We note that At'go^im = 1-47 x 10^^ [Hz] and At/i 4 o^m = 
0.831 X 10^^ [Hz] denote the band widt h for the AKARI WID E-S 
and WIDE-L band, respectively (see also lHirashita et alj2008l) . We 
derive 90pm and 140pm luminosity density from the observed flux 
density by multiplying 47rd|,/(l + 2 ) for each galax3[3. Here we do 
not consider fc-correction term as the AKARI/FIS band widths are 
wide enough, and its effect is negligible for our low-redshift galaxy 

sample (2 < 0.1). _ 

_ The IR-based SFRs (SFRir) are derived using iKennicuti 

( Il998h equation with the Kroupa IMF to be consistent with our 
UV- and Ha-based SFRs; i.e. SFRir=3.1 x 10"'‘'‘Lir [erg-s“i]. 
In Fig. [4] we compare SFRuv+ir(=SFRuv+SFRir) and dust- 
and aperture-corrected Ha-based SFRs. The scatter is not small 
(~0.3 dex), but these two completely independent measurements 
show a good agreement (with a median difference of ~0.01 dex), 

® We note that 1117 out of 1200 galaxies (~94%) within our AKARI- 
SDSS-GALEX sample have photometry at both 90pm and 140pm. The 
remaining 73 galaxies (~6%) are detected only at 90pm. We do not esti¬ 
mate Ljr for those detected only at 90pm, and they are not used in the 
following analyses. 
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Figure 5. The Ha extinction (Ajja) derived from Balmer decrement plotted 
against the observed SERrci/SFRuv ratio for all SF galaxies in our SDSS- 
GALEX sample. Dust extinction correction is not applied when comput¬ 
ing SFRhc« obs/SFRuv obs- The dashed line shows the best-fitted relation 
computed using the galaxies within the range of —0.5<log(Ho/UV)<1.5. 
There exists a positive correlation between the two quantities (as expected), 
but there remains a considerable scatter around the best-fitted relation. In 
the bottom-right comer of this plot, we also show the typical (median) error 
bars for individual data points derived from the flux uncertainties. 

considering the uncertainties associated with AKARI fluxes, as 
well as the Ha aperture/dust correction. We note that the catalogue 
constructed here will be published on the AKARI website. 

3 RESULTS 

3.1 Simple conversion from Ha/UV ratio to Ana 

Our first step is to check if there is really a positive correlation 
between the Ana and Ha/UV ratio, using all our star-forming 
galaxy sample. The Ha luminosities and FUV luminosity density 
are both proportional to recent SF activity, whereas these two indi¬ 
cators are affected by dust extinction at a different level. Because 
Ha line is less sensitive to dust extinction effects, galaxies with 
higher Ha/UV ratio are expected to be dustier. We plot in Fig. 
the Ana (derived from Balmer decrement) of all SF galaxies in 
our SDSS-GALEX sample against their SERHa.obs/SFRpuv,obs 
ratitfl As expected, there exists a positive correlation between Ana 
and SFRHc«,obs/SFRFuv,obs with the best-fitted relation of: 

Ahc = 0.743 X Iog(Ha/UV) -b 0.676, (5) 

where Ha/UV denotes the observed Ha-to-FUV SER ratio 
(SERHa,obs/SERpuv,obs). At the same time, it is also notable 
that there exists a substantial scatter around the best-fitted relation 
(ct~ 0.4 mag). Therefore, the data suggest that the dust extinction 
level of galaxies can (roughly) be approximated by Ha/UV ratio, 

® We note that all the analyses presented in this paper are based on the 
SFRna.obs/SFRpuv obs ratio for convenience, but our results are un¬ 
changed even if we instead use e.g. Lna obs/r^Tpuv obs (without con- 
veifing their luminosities to SFRs). 
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Figure 6. {Left): The Ahq: versus HqAJV plot (same as Fig. [5), colour-coded according to the average stellar mass at each position. We here applied 90x90 
gridding, and compute average log(M* /Mq) in each grid. We note that we request a minimum sample size of Afpix=4 for computing average stellar mass. 
We note that we apply this strategy when making similar plots in the remainder of the paper. (Right): The Ajja versus SFRjja/SFRuv diagram for different 
stellar mass bin as indicated in the plot, demonsteating that dust extinction level is strongly dependent on stellar mass of galaxies even at fixed Ho/UV ratio, 
in the sense that more massive galaxies are more highly obscured by dust. We apply 60x60 gridding (instead of 90x90) to create this plot, because of the 
smaller sample size in each sub-panel. The dotted line indicates the best-fitted relation for the total sample (Eq. 5). 


but it may be too uncertain to derive dust extinction by simply ap¬ 
plying the best-fitted zlna-HaAiV relation shown as the dashed 
line in Fig. In the remainder of this paper, we will examine the 
origin of this scatter, and attempt to reduce the systematic uncer¬ 
tainty in deriving the dust extinction from UV-optical information. 

In Fig. |5] a fraction of galaxies show log(Ha/UV)<0 (i.e. 
SFRhc«<SFRuv)- One may find it surprising, because it is ex¬ 
pected that SFRna.obs must be higher than SFRFuv,obs un¬ 
der the assumption that Ha is less sensitive to dust ex t inctio n 
(^FUv>^Ha)- However, as discussed by e.g. iLee et al] l l2009l) . 
UV-based SFR can exceed Ha-based SFR for various reasons, par¬ 
ticularly for low-mass galaxies. Indeed, most of the galaxies ex¬ 
hibiting log(HaAJV)<0 in our sample turned out to be low-mass 
galaxies (see Section 4.2). We emphasize that the aim of this paper 
is to empirically link the observed Ha/UV flux ratio to the Ha/H/3 
line ratio (hence dust extinction level), and so we do not discuss 
this issue further in detail in the remainder of this paper. 


3.2 Dependence on stellar mass 

Recent studies have shown that dust extinc tion of galax i es are 
strongly correlate d with stellar mass (e.g. iGam & Best! bOld: 
IReddv et al.ll2010n . Therefore we first investigate stellar mass de¬ 
pendence of the Ano-Ha/UV relation. In the left panel of Fig.| 6 ] 
we show how the average stellar mass of galaxies change on the 
Ano-Ha/UV scattered plot. Fig. is the same plot as Fig. 0 
but colour-coded according to the average stellar mass computed 
at each point (redder colours indicate higher M*). We here di¬ 
vide our sample into 90x90 grid and compute average stellar mass 
(log(M*/M 0 )) in each grid. We request a minimum sample size 
of Wpix = 4 for computing the average stellar mass, and pixels 
containing the smaller number of galaxies (i.e. with Wpix < 3) are 
not shown. 

A visual inspection of the left panel of Fig. reveals that 
more massive galaxies tend to be more highly obscured by dust 
at fixed Ha/UV ratio. By comparing the most massive galaxies 


(log(M*/M 0 ) ^ 10.5) and low-mass galaxies (log(M*/M 0 ) g 
9.0) in our sample, we find that the systematic difference in terms of 
Ana at fixed Ha/UV ratio can reach > 1.5-mag level at maximum. 
This result suggests that it might be misleading to blindly convert 
Ha/UV ratio into Ana for individual galaxies without consider¬ 
ing the stellar mass difference. In the right panel of Fig. we plot 
Ana against Ha/UV ratio by dividing the sample into five stellar 
mass bins. This plot also demonstrates that the location of galaxies 
on the Ahq!-H a/UV diagram is largely dependent on their stellar 
mass: high-mass and low-mass galaxies dominate completely dif¬ 
ferent regions on this diagram. 

It is also interesting to point out that the “slope” of the Ana- 
Ha/UV relation can change with stellar mass. The non-zero slope 
of the Ano-Ha/UV relation at fixed stellar mass seen in Fig. [^im¬ 
plies that the dust extinction cannot be determined solely from stel¬ 
lar mass alone, and this result motivates us to establish an empirical 
calibration for deriving more realistic estimates of Ana using the 
Ha/UV ratio and stellar mass of galaxies. 

We here assume the following conversion equation: 

Ah. = o(log M*) X log(Ha/UV) + 6 (log M*), ( 6 ) 

where a(log M^,) = ai log M* + 02 and b(log M*) = 61 log M* + 
62 - We fit all SDSS-GALEX sample with the above equation, and 
the resultant a(log M*) and &(log M*) are derived as follows: 

a(log M*) = 0.210 X log M* - 1.597, (7) 

6 (log M*) = 0.493 X log M* - 4.121. ( 8 ) 

The stellar mass trend reported here could partly be explained by 
a s tellar mass d e pende nce of the star formation history. As shown 
by IWuvts et al.l ( l2013l) . an intrinsic Lh./Tuv ratio evolves with 
galactic age and/or SF history of galaxies, in the sense that Ha/UV 
ratio rapidly declines with increasing age fo r those having exp o- 
nentially declining SF history (see fig. 3 in IWuvts et ahlEoiAn . 

^ IWuvts et all J2013l) show the evolution of the TH./T 28 OO ratio (not 
Th./Tfuv), bat the situation is qualitatively the same. Since contribu- 
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Figure 7. The same plot as Fig.[^ but colour-coded based on the average EWhck at each pixel (left panel). We also show how the Ahq;-Hq;/UV correlation 
changes with different EWhq (right panel). Galaxies with higher EWhck (i-e. with younger ages) tend to be less obscured by dust at fixed HaAJV ratio. 
Similar to Fig.[^ we apply 60x60 gridding to create the right-panel of this plot. 


Although this effect is expected to be small , as long as the SFR 
does not drop very rapidly with time (see e.g. iHao et al]|201lh . our 
result may suggest that the SF history of massive galaxies are more 
like exponentially declining with shorter time scale, resulting in a 
FUV flux excess (compared with Ha) due to the contribution from 
lower-mass stars. 

Another possibility is that the stellar mass trend seen in 
Fig. [6] can be explained by a stellar mass dependence of dust 
geometry within the galaxies. It is expected that Hq emission 
originates from star-forming regions (near O-type stars) with 
enhanced level of dust extinction due to the effects of opti¬ 
cally thick short-lived birth clouds. As a result, nebular emis¬ 
sion (including Ha) suffers extra extinction correction com¬ 
pared to_the_conhnuum_liglh^_at^_the_sarne_wavelength range 
(e. g. I Calzetti, K inn ey, & Storchi-BergmannI ri994h . As illustrated 
bv IPrice et alJ (20141), the more the SF regions are distributed uni¬ 
formly over the galaxies, the difference between stellar and neb¬ 
ular extinction level becomes smaller. On the contrary, if stars are 
formed in a compact region (e.g. in the central core of galaxies), Ha 
emission would more severely suffer from the dust extinction ef¬ 
fects than continuum light at the same wavelength. Our result sug¬ 
gests that more massive galaxies tend to have higher level of extra 
extinction towards the nebular region (Anebuiar>^cont), whilst SF 
regions in low-mass systems tend to be more widely spread over 
the galaxies (hence Anebuiar ~ Acont). We note that this point will 
be further discussed in Section 4.4. 

We also comment that we cannot completely rule out the pos¬ 
sibility that the observed trend is (partly) produced by the aperture 
correction effect—we recall that the Ha fluxes are measured with 
a limited size of aperture. Although we believe that the effect is 
small, this point will be discussed more in detail in Section 4.1. 


tion from low-mass stars to the FUV luminosities should be smaller than 
NUV, the age dependence of the ratio could be even smaller 

than that of iiHa/F/NUV- 


3.3 Dependence on EWhq 

We pointed out in Fig. that there remains a positive correlation 
between A-tia and Ha/UV at fixed stellar mass. This implies that 
we cannot fully determine the dust extinction of galaxies with stel¬ 
lar mass alone. In this section, we will investigate the Ana ver¬ 
sus Ha/UV relation as a function of another observable, EWhci, 
which is generally more directly linked to the age of galaxies (i.e. 
younger galaxies have higher EWna). We note that EWhc* indi¬ 
cates the rest-frame Ha equivalent widths in the remainder of this 
paper. 

The left panel of Eig. [7] shows the Ana-Ha/UV diagram, 
colour-coded based on the average EWhc« computed at each point 
(i.e. redder colours indicate higher EWna). As we performed in 
Fig.H we apply 90x90 gridding and compute average EWh^ in 
each pixel, by requesting a minimum sample size of A^pix = 4 in 
each pixel. This plot demonstrates that galaxies with higher EWhc* 
tend to show lower dust extinction at fixed Ha/UV ratio (or equiv¬ 
alently, galaxies with higher EWhq tend to show higher Ha/UV 
ratio at fixed Ahq), which can probably be explained by the same 
reason as discussed in the previous section: i.e. galaxies with higher 
EWho are expected to have younger stellar age with little contribu¬ 
tion from old stars to the EUV luminosities. In fact, galaxies with 
EW>50A runs through the (0, 0) point on Fig. [7] whilst low-EW 
galaxies show significant offsets towards higher Ahc- This result 
again suggests that a simple conversion from Ha/UV to Ahc could 
be too simplistic. 

Following the procedure that we adopted when deriving the 
stellar mass dependence of the Ana-Ha/UV relation in Sec¬ 
tion 3.2, we here attempt to obtain an empirical calibration for de¬ 
riving Ahc« from EWhc« and Ha/UV ratio. We assume a similar 
form of conversion equation as we did in Section 3.2: 

Ahc = a(log EWhc) X log(Ha/UV) -F 6(log EWhc), (9) 

where a(logEWHa) = ai logEWno -F 02 and fe(logEWHa) = 
bi log EWho -Fb 2 . We fit all SDSS-GALEX sample with the above 
equation, and the resultant a(logEWHc«) and 6(logEWHa) are 
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Figure 8. (Left): Dust extinction (Ajja) derived from Balmer decrement plotted against the rest-frame EWhh- The colour coding indicates the stellar mass at 
each pixel, as we performed in the previous plots. (Right): The same plot for different stellar mass bin as indicated in the plot. It is evident that the location of 
the Aho versus EWhq relation is strongly dependent on stellar mass. Interestingly, the slope of the Ajja versus EWhq relation changes with stellar mass: 
Ajjo increases with increasing EWjjq for massive galaxies, whilst the relation is nearly flat for low-mass galaxies. 


derived as follows: 

a(logEWHa) = 0.101 X log EWh£« +0.872, (10) 

6 (logEWHa) = —0.776 X logEWna + 1.688. (11) 

This conversion equation can be an useful tool for predicting dust 
extinction, particularly when stellar mass estimate is not available. 
Another advantage of this prescription is that it is described only 
with observed quantities (EWhc« and Hq/UV without correction), 
and therefore it does not suffer from uncertainties regarding stellar 
mass estimates accompanied by SED fitting. 

The stellar mass dependence of the Ana-Ho/UV ratio pre¬ 
sented in Section 3.2 is based on all SF galaxies selected from 
SDSS sample using BPT diagram (Section 2.1). Flowever, we want 
to stress that the definition of SF galaxies can be different from 
studies to studies: e.g. the sample is usually biased to higher EWhc* 
galaxies in the case of high-^ studies (hence biased to higher 
specific-SFR galaxies with young stellar population). We will dis¬ 
cuss later in Section 4.5 whether our new prescription can be appli¬ 
cable to high-redshift galaxies. 


3.4 Predicting dust extinction without UV information 

The main aim of this paper is to establish an empirical link be¬ 
tween the Flof/UV luminosity ratio and dust extinction properties 
(as presented in Section 3.1-3.3), but our analyses presented in the 
previous sections demonstrate that Ana is clearly dependent on 
stellar mass and EWhq. In Figs. [6] and [7] it can be seen that the 
“fixed-mass line” and “fixed-EW line” show different slope. This 
motivates us to construct another prescription to predict Aho for a 
given mass and EWhq. 

In the left panel of Fig. [ 8 ] we plot Aua against EWhq. The 
colour-coding indicates average stellar mass at each pixel (redder 
colours indicate higher stellar mass). There seems to be a general 
trend that Ahq declines with EWna, but obviously, there is a large 
variation in the dust extinction properties at a given EWhci- It is 
clear from Fig. [ 8 ] that there is an universal trend that more mas¬ 
sive galaxies tend to be more highly obscured by dust at fixed 
EWhq. In addition, it is interesting to note that, at fixed stellar 


mass, Ahc« sharply increases with increasing EWna for massive 
galaxies, whilst this trend is not visible for low-mass galaxies (Ana 
always stay at gO.5 mag regardless of EWhq: see the right panel 
of Fig. [§}. This result suggests that the nature of SF galaxies is 
strongly dependent on stellar mass. 

We can now establish an empirical prescription for deriving 
Ahci from M* and EWna- We assume a similar form of conversion 
equation as we did in the previous sections (but this time we do not 
require UV information): 

Ah. = a(logM*) x log(EWH.) + b(logM*), (12) 

where a(log M*) = ai log M*+a 2 and 6 (log M^,) = bi log M* + 
62 - We fit all SDSS-GALEX sample with the above equation, and 
the resultant a(log M*) and fo(log M*) are described as follows: 

a(log M*) = 0.096 x log M* - 0.717, (13) 

&(log M*) = 0.538 X log M* - 4.745. (14) 

Our results suggest that high-mass and low-mass galaxies 
show different behaviour on the Ah. versus EWh. diagram, prob¬ 
ably indicating that their dust extinction properties are different. It 
is possible that massive galaxies tend to become more like nuclear 
starbursts with increasing EWh. (hence dusty starbursts with high 
specific-SFR), whilst low-mass galaxies tend to form stars over the 
galaxy discs, regardless of their specific-SFR. 

3.5 Summary: a new recipe for predicting dust extinction 
correction in the absence of Balmer decrement 

We have demonstrated that, as expected. Ah. derived from 
Balmer decrement shows positive correlation with the observed 
SFRh./SFRuv ratio (Section 3.1). This is probably a natural con¬ 
sequence reflecting the fact that Ha is less sensitive to dust extinc¬ 
tion than UV continuum light, consistent with what we had initially 
expected. However, we pointed out that there remains a large scatter 
around the Ah.-HqAJV correlation, and therefore we suggest that 
predicting dust extinction of individual galaxies using only HaAJV 
ratio could be highly uncertain. 

In this paper, we have shown that the scatter around the Ah.- 
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Table 1. Summary of the new recipes for predicting dust extinction proposed here and the conesponding sections in the paper. 


Recipe ai 02 bi 62 crnns [mag] Section 


(1) Ha/UV only. — — 0.743±0.004 0.676±0.002 0.386 §3.1 

(2) mass + Ha/UV.. 0.210±0.006 -1.597±0.068 0.493±0.003 -4.121±0.026 0.282 §3.2 

(3) EW + Ha/UV.. 0.101±0.014 0.872±0.020 -0.776±0.006 1.688±0.008 0.338 §3.3 

(4) mass + EW.. 0.096±0.007 -0.717±0.073 0.538±0.011 -4.745±0.104 0.317 §3.4 



-3 


g 10 11 

log(U*/HS) 



3 10 11 


log(M*/MO) 



-3 


g 10 11 

log(M*/MO) 



10 

log(U*/MO) 


(d) Mass 


EW(Hd)J 


Figure 9. Comparison of the new recipes for predicting Aho proposed in 
this paper. Here we show the difference between the Ajja derived from 
Balmer decrement and those derived using the new recipes. The top panel 
shows the result for the case using only Ha/UV ratio. The scatter becomes 
smaller in all cases (panels [b]-[d]), although they are not necessarily dra¬ 
matic improvements (see also Table 1). We also show the FIR-detected 
gala xies (AKARI sources) wi th red circles, as well as 2 : ~ 2.5 galaxies 
from IShimakawa et alj i2015bll with blue circles. 


Ha/UV correlation can be reduced by combining stellar mass 
and/or EWhq, and we proposed new recipes for predicting dust 
extinction correction using a combination of Ha/UV, stellar mass, 
and EWhci- The first approach is to use stellar mass and Ha/UV ra¬ 
tio (Section 3.2). The second approach is to use EWh^ and Ha/UV 
ratio (Section 3.3), which does not require stellar mass estimate, 
so that it is particularly useful when multi-band photometry is not 
available. The final approach presented in this paper is to use stel¬ 
lar mass and EWhc« (Section 3.4). This approach has an advantage 
in the absence of UV information. In Table 1, we summarize the 
proposed recipes outlined here. 

Overall, which is the best way to predict dust extinction among 
the proposed recipes? From Table 1, it can be seen that we can re¬ 
duce the scatter around the best-fit relation in all cases (by ~15- 
30%), with the “mass - 1 - Ha/UV” approach being the best. Al¬ 
though they are not necessarily dramatic improvements, we stress 


that these new recipes would be able to provide more realistic dust 
extinction levels for individual galaxies than the simple Ha/UV 
approach. In Fig. we show the difference between Ahc from 
Balmer decrement and those from our new recipes, as a function of 
stellar mass. In the case of simple Ha/UV method (panel-[a]), there 
clearly remains a systematic uncertainty: we tend to underestimate 
A-tia for massive galaxies, while we tend to overestimate Aua for 
low-mass galaxies. Interestingly, for low-mass galaxies, a constant 
~0.5-1.0 mag correction would be more realistic than applying the 
Ha/UV approach (see panels for low-mass galaxies in Fig. B- 
In Fig. we also plot our AKARI FIR-detected galaxies as 
well (red symbols). Their dust extinction level seems to be signifi¬ 
cantly underestimated by the simple Ha/UV approach (as they tend 
to be massive galaxies), but the situation is improved with the new 
recipes. An important message of this paper is that we can predict 
a reasonable dust extinction correction (for Ha) using UV-optical 
information alone, even in the absence of H/3 lines or deep FIR 
photometry. 


4 DISCUSSION 

4.1 Systematic effects regarding aperture correction 

The most important results presented in this study is that the A^a- 
Ha/UV relation is strongly dependent on M* and EWhq. A ma¬ 
jor concern when interpreting our results is the impact of aperture 
correction. As mentioned in Section 2, Ha flux measurements are 
made with a limited size of SDSS fibre (3” in diameter). We recall 
that the fibre can only trace 10-20% of the total light (at r-band) 
in extreme cases (see Fig. 13, and so one may be worried that there 
remains a large uncertainty in the total Ha flux. 

In this study, we applied aperture correction to derive Ha to¬ 
tal luminosity based on the difference between fibre- and total- 
magnitudes measured at r-band continuum (see Section 2.1). If 
massive galaxies tend to have quiescent bulge component in their 
central part, and if those galaxies tend to have active star formation 
in their outskirts, we could underestimate Ha aperture correction: 
i.e. SFRhq could be underestimated in this case. However, we be¬ 
lieve it is not likely the case, because our results do not change even 
if we use SDSS tt-band photometry (which is expected to reflect 
more recent SF activity than r-band) to estimate aperture correc¬ 
tion. 

In Fig.[T3 we plot the difference between the aperture correc¬ 
tion values measured with u-band and r-band (Au — Ar) as func¬ 
tions of various galaxy properties. The stellar mass dependence is 
only -^0.2-mag level at most, showing that the strong stellar mass 
dependence reported in Fig. [6] (Section 3.2) cannot be fully ex¬ 
plained by the aperture correction effect. In Fig. [TO] it can be seen 
that the trend is most significant for EWhc*: the Au — Ar value 
declines with increasing EWna- This would indicate that galax¬ 
ies with higher EWhc* tend to be forming stars in more compact 
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Figure 10. Difference of aperture correction derived from u-band and r-band (Ait — Ar) plotted against stellar mass, SFR, EWho. and Ajja (from left to 
right). This plot demonstrates that the aperture correction factors derived from u-band and r-band agree reasonably well, and therefore we expect that the 
aperture correction issue would not strongly bias the results. We note that the trend is most significant for EWho. indicating that u-band light tend to be more 
strongly concentrated in the central region for higher-EW galaxies. Nevertheless, we stress that our conclusion does not change even if we use the aperture 
correction derived from it-band data. 


regions (hence ii-band light is more centrally concentrated), but 
the difference is still too small to fully explain the trend we have 
shown in Fig.|7](Section 3.3). We note, however, that u-band is not 
necessarily a perfect way to trace very recent star forming activ¬ 
ity (compared with Ha), and so we cannot completely rule out the 
possibility that even u-band light distribution might differ from that 
of Ha emission. Furthermore, we need to assume in this paper that 
EWho and Ana derived from fibre spectroscopy can be applica¬ 
ble to outer part of galaxies. This is inevitable as far as we rely on 
the SDSS data, and it is essential to make observations to map Ha 
and H/3 emission over the galaxies with large aperture, in order to 
completely verify our results. 

We recall that the adopted aperture correction value is corre¬ 
lated with EWho as reported in Eig.j^ i.e. we need to apply larger 
aperture correction for lower YANna galaxies. We therefore caution 
that the levels of aperture correction also change on the Ana versus 
Ha/UV plot, in the sense that galaxies located at the upper enve¬ 
lope of the ylna-Ha/UV relation tend to require larger aperture 
correction (Eig. [ID- Nevertheless, as discussed in Section 2.1, we 
attribute this trend to the physical reason that galaxies with higher 
EWho tend to be more compact, and we believe that the trend 
between EWna and the levels of aperture correction reported in 
Eig-HUdoes not make any artificial bias on our conclusion. 


4.2 Can we predict Ana with stellar mass alone? 


Recent studies suggest that stellar mass is the most important pa¬ 
rameter which determines the dust extinction levels for SE galax¬ 
ies (e.g. lGam & Bestiboiol) . Using our own datasets, we also con¬ 
firm that the scatter around the best-fit ^Ha-stellar mass rela¬ 
tion is fairly small (0.32 mag): this is actually even smaller than 
the Ha/UV approach (~0.38 mag), and is comparable to those 
achieved with our new recipe (see Table 1). Eurthermore, the cor¬ 
relation between dust extinction an d stellar mass is reported to be 
unchanged at least out to z ^ 1.5 l Gam et aDl2010l : ISobral et al.l 


l2012tlibar et al.ll2013l : lDommguez et alJl2013*1) 7 and so stellar mass 

can be an useful indicator for a rough estimate of dust extinction 
of galaxies at all redshifts in the absence of the B aimer decrement 
or deep EIR photometry (but see IOtedl201^ : IPannella et al.ll201^ : 
IShimakawa et aklboiSah . However, as demonstrated in this work, 
it might be too simplistic to blindly apply the M*-dependent dust 
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Figure 11. The Ana versus Ha/UV plot (same as Figs. 5-7), with the grey¬ 
scale indicating the average aperture con'ection value for galaxies at each 
pixel. Galaxies located top-left corner of this plot tend to require larger 
aperture correction. We recall that galaxies located at the top left coiuer 
of this plot tend to have lower EWho (see Fig. 7). Therefore, the trend 
seen in this plot is actually equivalent to Figs. 2 and 7, where we showed a 
decreasing trend of aperture correction with EWna- 


extinction correction. In this section, we discuss how well we can 
really predict dust extinction from stellar mass alone. 

In the left- and middle-panel of Eig.[T3 we plot EWna and 
SER as a function of stellar mass. The colour-coding indicates the 
average Ana derived from B aimer decrement at each point (red¬ 
der colours indicate higher dust extinction). The left panel is es¬ 
sentially equivalent to Fig. where we showed the stellar mass 
dependence of the Hna-EWHa relation. This plot further demon¬ 
strates that more massive galaxies tend to be dustier on ave rage 
(consistent with many other studies; e.g. lWhitaker et all2012h . but 
it is also important to note that, at fixed M*, galaxies with higher 
EWho tend to be more highly obscured by dust. On the other hand, 
the middle panel of Fig. [7^ shows how the average dust extinction 
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Figure 12. (Left): The EWho versus stellar mass diagram, colour coded based on the dust extinction (Ajja)- The fact that a fixed-AHa line does not run 
completely perpendicularly on this diagram suggests that we can provide a more reliable estimate of dust extinction by using stellar mass and EWua ■ (Middle): 
A similar plot to the left panel, but here we show the dust extinction levels of galaxies on the SFR-M* diagram. This plot again shows that more massive 
galaxies tend to be dustier, but the behaviour is slightly more complex compared to the left-hand panel, implying that EWuct can be a better tool for estimating 
dust extinction of galaxies. (Right): Dust extinction versus stellar mass plot, colour coded based on the observed SFRhc«/SFRuv ratio. This plot further 
demonstrates that predicting dust extinction from stellar mass alone can be too simplistic, whilst we can reduce the systematic uncertainty by using stellar 
mass and FI«/UV ratio. This is the most important finding of this work. 


of galaxies changes on the SFR-M* diagram. Again, more massive 
galaxies tend to be dustier, whilst the behaviour of a fixed-A hq line 
on the SFR-M* diagram is rather more complex compared with the 
EW-M* plot (left panel). It is interesting to note that a more simple 
observable (EWna) can be a better tool for predicting dust extinc¬ 
tion of galaxies—and this is why we exploit EWhc* (a proxy for 
the age of stellar population) rather than SFRs when we establish 
an empirical calibration in Section 3. 

Finally, in the right panel of Fig. d we show Ana against 
stellar mass. Here we apply colour-coding based on the Ha/UV ra¬ 
tio as we investigated in this paper. This plot is equivalent to Fig.[^ 
but more directly demonstrates our important result that a simple 
stellar-mass dependent dust extinction correction is not a perfect 
way: by adding Ha/UV ratio we can have a more realistic dust ex¬ 
tinction correction in the absence of Balmer decrement and/or deep 
FIR photometry. 


4.3 Metallicity dependence 

Dust is formed by metals, and so the du st extinction may be 
correlated with metallicity of galaxies (e.g. IXiaoetalJ[20l% . In 
Fig.d we plot Ha dust extinction derived by Balmer decrement 
for our SDSS-GALEX star-f orming galaxies again st their (gas- 
phase) metallicity derived by iTremonti et al.l l l2004) . The colour 
coding of this plot indicates average stellar mass at each point (red¬ 
der colours indicate higher M*). 

It is clear from Fig.ll3lthat there is an overall trend that dust 
extinction level increases with increasing metallicity. It is also evi¬ 
dent that both dust extinction and metallicity increase with increas¬ 
ing stellar mass. At fixed metallicity, dust extinction level slightly 
increases with stellar mass. Similarly, at fixed Ahc, metallicity 
also increases with stellar mass. At fixed stellar mass, on the other 
hand, it looks as if dust extinction level sharply decreases with in¬ 
creasing metallicity. This may sound surprising, because the trend 
is completely opposite to the overall dust-metallicity correlation. 
We note, however, that this is just an average behaviour—we find 
that individual data points at fixed stellar mass do not align so 
tightly on the fixed-mass line in Fig. Neveretheless, it is in- 



Figure 13. Dust extinction versus gas-phase metallicity for our SDSS- 
GALEX star-forming galaxy sample, colour-coded based on the average 
stellar mass at each point. We applied 60 X 60 gridding for making this plot. 
There is a general trend that the dust extinction level increases with in¬ 
creasing metallicity. It is also notable that the dust extinction increases with 
stellar mass at fixed metallicity, and at the same time, metallicity increases 
with stellar mass at fixed Ajja. 

teresting to point out that a small variation (scatter) around the 
mass-metallicity relation is related to the variation of dust prop¬ 
erties of galaxies at fixed stellar m ass (see consistent results for 
high -2 galaxies bv IZahid et al.ll2014l) . 

Detailed modeling of the chemical evolution and dust forma¬ 
tion/destruction processes within galaxies is beyond the scope of 
this paper, but we comment that our results are basically unchanged 
even if we replace the stellar mass parameter with the metallicity. 
However, we recall that the main aim of this paper is to establish an 
empirical link between dust extinction properties of galaxies and 
some quantities which can easily be obtained. Deriving metallicity 
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Figure 14. Compaiison between the Ha (nebular) dust extinction from Balmer decrement and UV (continuum) dust extinction measured with IRAJV ratio, 
for AKARl FIS detected galaxies. The colour coding of the data points represent the stellar mass (left panel) and the rest-frame EWhc« (right panel): i.e. 
redder colours indicate higher M* and higher EWho in the left and right panel, respectively. The three dotted lines correspond to the extra extinction factor of 
/ = 0.44, 0.7, and 1.0 (in the case of Calzetti extinction law). This plot demonstrates that more massive galaxies and/or lower EWhq galaxies tend to have 
higher extra extinction towards the nebular regions (i.e. smaller / value). 


of galaxies is usually a more complex issue (and it requires detec¬ 
tion of several faint emission lines), hence we decided to use stellar 
mass (rather than metallicity) when we derive the empirical cali- 
hration in Section 3. 


4.4 Stellar mass and EWhq dependence of the extra 
extinction towards nebular regions 

It is suggested that there is a need for an extra extinction correction 
to nebular emission lines (including Ha) compared with contin¬ 
uum light at the same wavelength. This is because nebular emis¬ 
sion lines are originated in the dust-rich birth cloud, whilst other 


(see e.g. Calzettil 19971: 

Chariot & Fall! 

2000 

: ICid Fernandes et al.l 

l2005l:IWuvtsetal.ll2011 

: Steidel et alj 

2014 

). Using a sample of 

local star-forming galaxies. ICalzetti et alj 
colour excess of the stellar continuum, Es{ 
that of nebular emission lines, Sgas {B — V 

20001) argue that the 
B — V), is linked to 
, with Es{B — U) = 


0.44 X E^as{B — V) (i.e. extra extinction factor, / = 0.44). In 
fact, the amount o f this extra extinctio n factor (/) is still under de¬ 
bate. For instance. I Wuvts et'H] ( 1201 3l) show a polynomial function 
to derive Hextra from the continuum extinction (Acont), hence the 
/ value may not necessarily be a constant value for all types of 
galaxies. It is also suggested by recent studies that the typical / 
value m ay be higher (/ = 0.5 — 1.0) for high-redshift galaxies 


(see e.g. Erbet^lj 20061 : iReddv et'd] 20101 : iKashino et al.ll2013l : 
I Valentino et al.ll201^ 

Taking advantage of the AKARl FIR photometry (which cap¬ 
tures the peak of dust emission), we compare in Fig. [14] the dust 
extinction derived from Balmer decrement (Ahc«; i.e. nebular ex¬ 
tinction) and the IR/UV ratio (Auv; i.e. continuum extinction). It 
can be seen that there is broadly a positive correlation between the 
two independent measurements of dust extinction properties, but 
the left panel of Fig.[T4|also demonstrates that more massive galax¬ 


ies tend to have higher Ahq than low-mass galaxies at fixed Auv- 
Similarly, galaxies with lower EWh^ tend to have higher Ahq at 
fixed Auv (see right panel of Fig.[T4t. 

In Fig. [H we show three dotted lines corresponding to 
/ =0.44, 0.7, 1.0, after properly converting Auv to the continuum 
extinction at Ha. We find that most of our galaxy sample (AKARI- 
detected star-forming galaxies) are distributed in the range between 
/=0.44 and /=1.0, and at the same time, we also find an interest¬ 
ing trend that more massive galaxies tend to have smaller / value: 
i.e. higher extinction towards the n ebular regions. Interestingly, a 
very recent study by IPuglisi et al.l l l2015l) made a similar sugges¬ 
tion for high-redshift galaxies (see their fig. 10). We therefore sug¬ 
gest that the stellar mass (or EWhq) trend of / value reported in 
our current study may hold over the cosmic time. Assuming that 
the difference in / va lue is due to dust geometry within the galax¬ 
ies l lPrice et alj|2014h . our results suggest that high -2 SE galaxies, 
which are often reported to have / ~ 1, are more like low-mass 
galaxies in the low-redshift Universe, which tend to be forming 
stars over the disc. 

We caution that the stellar mass and/or EWhq dependence 
of the / value reported here is implicitly incorporated when es¬ 
tablishing our empirical recipes for predicting dust extinction lev¬ 
els of galaxies in Section 3. For example, for massive galaxies, 
our Ha/UV methods automatically assume small f values. This 
might be misleading, particularly if high-redshift galaxies tend to 
have larger f value comparable to low-mass galaxies in the local 
Universe—in this case we might over-estimate the dust extinction 
levels of (massive) high -2 galaxies if we simply applied our lo¬ 
cal calibration. However, unfortunately, high -2 galaxy sample cur¬ 
rently available is too small to test this possibility (see blue symbols 
in Fig.|9ll. Also, by considering the situation that the typical ”/” 
value reported for high -2 galaxies significantly vary from studies 
to studies, it would be impossible at this moment to conclude the 
above possibility (see also Section 4.5). 
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We note that we have used ICalzetti et alj l l2000h extinction law 
to convert Auv to the (continuum) extinction level at Ha wave¬ 
length, to be consistent with all the analyses presented in this pa¬ 
per. If we instead used other extinction curves, then the absolute 
/ values would change by a factor of ~2x at maximum: e.g. 
more specifically, the slope of the dotted lines on Fig. M would 
be reduced by ~10% when we apply the Milky -way dust extinc¬ 
tion curve of ICardelli. Clayton, & Mathij ( Il9^ ), and by -^50% i n 
the case of the SMC-type extinction curve of loordon et alj l l2003h . 
Nevertheless, the relative difference in / values between high-mass 
and low-mass galaxies does not change as long as we assume the 
same type of extinction curve for all galaxies in the sample. We 
note that this assumption itself is uncertain, and the shape of ex¬ 
tinction curves may depend on g alaxy properties such as (specific) 
SFRs (see e.g. (Wild et alJl201 ih . However, unfortunately, it is im¬ 
possible to assess the shape of extinction curve (and its stellar mass 
dependence) with the current dataset alone. 


4.5 Application to high-redshift galaxies 

The main aim of this study is to construct an empirical calibration 
to predict dust extinction from the observed Ha/UV ratio (com¬ 
bined with some other galaxy properties). The new recipes de¬ 
veloped in this paper should be an useful tool particularly when 
B aimer decrement or deep FIR photometry is not available, which 
is usually the case for high-redshift studies. In this final subsec¬ 
tion, we discuss if our new recipes can really be applicable to high- 
redshift galaxies. 

Although it is becoming easier to have good quality NIR spec¬ 
tra for high-z galaxies, it is still very challenging to detect H/3 
lines for individual galaxies even with 8-m class telescope . We 
use 2 = 2.5 galaxy sample from IShimakawa et alJ l l2015bh . for 
which Ha and H/3 fluxes are both directly measured with deep 
spectroscopy with MOSFIRE on Keck. In Fig.|3 we plot our high- 
2 galaxy sample with the blue circles on top of the distribution of 
our local galaxy sample. Obviously, our high -2 sample currently 
available is by far too small to tell if our procedure works well 
for high -2 galaxies. Therefore, we stress that the applicability of 
our new recipe to high-redshift galaxies needs to be tested using 
larger high -2 galaxy sample with high S/N spectra in the future. 
A growing number of extensive NIR spectroscopic observational 
campaign with new facil ities are now under way (e.g. lSteidel et al.l 
I 2 OI 4 I : iReddv et^l2015h . and also, future thirty-meter class tele¬ 
scopes will of course deliver high-quality spectra for high-redshift 
galaxies. We believe that those new observations will revolutionise 
the situation and make it easier to access Balmer decrement for 
high -2 galaxies, allowing us to test if our new empirical recipes for 
predicting dust extinction properties established in the local Uni¬ 
verse really work for high-redshift galaxies. 


as EWho of galaxies. At fixed Ha/UV ratio, galaxies with higher 
stellar mass or lower EWho tend to have higher Ana, which could 
be explained by the stellar mass dependence (and/or EWna depen¬ 
dence) of the SE history: an intrinsic Ha/UV luminosity ratio is 
expected to rapidly decline with galac tic age, particularly when as¬ 
suming exponentially declining SER dWuvts et al.ll2013h . Another 
possibility is that more massive galaxies tend to have higher ex¬ 
tra extinction towards nebular regions, so that Ha is more heavily 
attenuated with respect to UV continuum light (see below). 

At fixed stellar mass. Ana is still positively correlated with 
Ha/UV ratio and with EWna, suggesting that predicting dust ex¬ 
tinction from stellar mass alone could be too simplistic. We have 
established an empirical calibration for predicting Ana using (1) 
stellar mass - 1 - Ha/UV ratio, (2) EWho + Ha/UV ratio, and (3) stel¬ 
lar mass -I- EWhu. By comparing the Ana value derived by Balmer 
decrement and those from the newly proposed methods, we find 
that our new methods work reasonably well (and successfully re¬ 
duce the scatter around the best-fitted relation by ~ 15-30%), with 
the “mass - 1 - Ha/UV” method being the best. Nevertheless, the other 
methods can also be useful in any future studies, particularly when 
stellar mass estimates and/or UV continuum information are not 
available. 

We find that Ana for high-mass galaxies increases with in¬ 
creasing EWho, whilst dust extinction does not increase with 
EWho for low-mass systems, suggesting a different nature of SE 
activity for galaxies with different mass. We use AKARI EIR 
data to derive UV extinction (Auv) using IR/UV ratio, and test 
the stellar mass and EWho dependence of the extra extinction 
correction factor towards the nebular region (/ value, defined as 
Esta.r{B — V) = / X i?gas(i3 — U)). We find an interesting hint that 
the / value is dependent on stellar mass or EWho, with more mas¬ 
sive galaxies or low EWhq galaxies having higher extra extinction 
towards nebular regions (i.e. smaller / value). Considering recent 
studies claiming higher /-value for high -2 galaxies, we argue that 
the dust geometry within high -2 SE galaxies resemble more like 
low-mass galaxies in the nearby Universe. 

An important caveat is that our empirical recipes for predict¬ 
ing dust extinction from Ha/UV ratio implicitly incorporate this 
stellar mass and/or EWho dependence of the / value: e.g. for mas¬ 
sive galaxies our recipes automatically assume small f value. If 
high-redshift galaxies tend to have larger f value than local coun¬ 
terparts with the same mass as often reported by recent studies, 
our empirical recipes established with local galaxy sample might 
over-estimate the dust extinction levels for high-redshift galaxies. 
However, the current high -2 galaxy sample (with measured Ha/H/3 
ratio) is too small to test if our new recipes can be applicable to 
high -2 galaxies. 


5 SUMMARY 

We presented an empirical calibration between the observed Ha- 
to-EUV ratio and dust extinction (Ana derived from Balmer decre¬ 
ment) for SE galaxies at 0.02 < 2 < 0.1 (selected with BPT dia¬ 
gram) using our SDSS(DR7)-GALEX(GR5) matched galaxy sam¬ 
ple, helped by AKARI EIR all-sky data. We confirmed that Ana 
increases with increasing Ha/UV ratio (as expected), but there ex¬ 
ists a considerable scatter around the Ano-Ha/UV relation. 

An important finding of this study is that the scatter around the 
Ano-Ha/UV relation is largely dependent on stellar mass, as well 
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